Using acoustic waves to induce high-frequency current oscillations in superlattices 
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We show that GHz acoustic waves in semiconductor superlattices can induce THz electron dy- 
namics that depend critically on the wave amplitude. Below a threshold amplitude, the acoustic 
wave drags electrons through the superlattice with a peak drift velocity overshooting that produced 
by a static electric field. In this regime, single electrons perform drifting orbits with THz frequency 
components. When the wave amplitude exceeds the critical threshold, an abrupt onset of Bloch-like 
oscillations causes negative differential velocity. The acoustic wave also affects the collective behav- 
ior of the electrons by causing the formation of localised electron accumulation and depletion regions, 
which propagate through the superlattice, thereby producing self-sustained current oscillations even 
for very small wave amplitudes. We show that the underlying single-electron dynamics, in particular 
the transition between the acoustic wave dragging and Bloch oscillation regimes, strongly influence 
the spatial distribution of the electrons and the form of the current oscillations. In particular, the 
amplitude of the current oscillations depends non-monotonically on the strength of the acoustic 
wave, reflecting the variation of the single-electron drift velocity. 

PACS numbers: 73.21. Cd, 73.50.Fq, 73.50.Rb, 73.23.-b 



I. INTRODUCTION 

Electrons in semiconductor superlattices (SLs) exhibit 
a wide range of nonlinear effects that are of fundamen- 
tal scientific interest and useful for applications in ultra- 
fast electronics^—. Many of these effects originate from 
the SL minibands, which enable electrons to perform 
THz frequency Bloch oscillations when a sufficiently high 
static electric field is applied along the SL axisi^&i 4 -. 
Bloch oscillations cause the electron drift velocity to de- 
crease with increasing electric field, which can trigger 
charge-domain oscillations accompanied by the emission 
of electromagnetic radiatio n 11 ' 13 . 

The frequency response of SL oscillators and de- 
tectors is limited by scattering processes including 
electron-phonon interaction s 7 ' 13 ' 14 . Surprisingly, though, 
phonons can serve as a powerful tool for enhanc- 
ing the electronic and optical properties of solid state 
devices^ 4 -^. For example, in 'SASER' SLs2£, analogous 
to the laser, the amplification of coherent sound waves 
now opens the way to acoustic control of miniband elec- 
tron transport. In related work, we recently demon- 
strated that a train of acoustic strain pulses can induce 
current in a SL by dragging miniband electrons through 
the device^ 7 -. 

In this paper, we show that a continuous GHz acous- 
tic wave can create complex THz electron dynamics 
in SLs, thus producing high-frequency current oscilla- 
tions even when no static electric field is applied. In 
a single-electron picture, there are two distinct dynam- 
ical regimes, depending on whether the energy ampli- 
tude, U, of the acoustic wave is greater, or less, than a 
critical value, U c , which depends on the SL parameters. 
For U < U c , the acoustic wave drags electrons through 
the SL, producing, in the presence of electron scatter- 
ing, a drift velocity, Vd, far higher than the speed of the 



wave itself. In this regime, the electrons perform peri- 
odic orbits in the rest frame of the acoustic wave. The 
orbital frequencies are more than an order of magnitude 
higher than the sound wave frequency. When U > U Cl 
the acoustic wave no longer drags electrons through the 
SL. Instead, there is an abrupt onset of Bloch-oscillation 
bursts, which makes Vd decrease extremely rapidly with 
increasing U . 

Our results demonstrate that miniband electrons 
driven by a GHz acoustic wave can attain a higher max- 
imal Vd, and have a larger negative differential veloc- 
ity (NDV), than when they are accelerated by a static 
electric field. Since high-frequency SL oscillators require 
high values of both parameter s 13 ' 14 , acoustic wave driv- 
ing could strongly enhance the performance of such de- 
vices. To investigate how the sudden onset of single- 
particle Bloch oscillations affects the collective dynamics 
of the electrons, we solve the Poisson and (drift-diffusion) 
current-continuity equations self-consistently throughout 
the SL. We use our results to determine the current 
through the SL as a function of U and time. Remarkably, 
the current exhibits self-sustained high-frequency oscilla- 
tions for all U. For U < U c , these oscillations originate 
from space-charge density waves, created and dragged 
through the SL by the acoustic wave. In this regime, as U 
increases the space-charge modulation strengthens, and 
so the amplitude of the current oscillations also increases. 
But when U increases above U c , the onset of Bloch os- 
cillations localises the electrons, thereby weakening the 
current oscillations. This electron localisation also leads 
to the formation of propagating charge domains, which 
coexist with the charge density waves. 

The complex single-particle and collective electron dy- 
namics induced by an acoustic wave demonstrate the po- 
tential of SLs for interfacing high-frequency electronics 
with the emerging field of phononic o 28 ' 29 , which is at- 
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FIG. 1: Solid [dashed] curve: Vd versus U (lower scale) or e 
(upper scale) calculated for a miniband electron driven by an 
acoustic wave only [or accelerated by a constant electric field, 
ksll/e, only]. Dotted line [arrow] marks U = U c [Uet]- Inset: 
schematic diagram of the SL layers and co-ordinate axes. 



tracting considerable interest in applied physics and en- 
gineering. In particular, our results indicate that SLs 
can be used to both detect and up-convert an applied 
acoustic signal. Generic features of the energy band 
transport processes created by a propagating wave po- 
tential also suggest ways to control transport through 
other spatially-periodic systems, including cold atoms in 
optical Iattices2£i2i. 

The structure of the paper is as follows. In Section 
HT1 we introduce the semiclassical equations of motion 
for a miniband electron in the presence of an acoustic 
wave and show that the electron dynamics and drift ve- 
locity depend critically on the wave amplitude. In Sec- 
tion |TTT1 we consider the collective electron dynamics and 
show that the acoustic wave triggers linear charge den- 
sity waves, nonlinear charge domains (when U > U c ), 
and self-sustained current oscillations even for very small 
wave amplitudes. Finally, in Section IIV1 we summarise 
our results and draw conclusions. 



II. MODEL OF SINGLE ELECTRON 
DYNAMICS 

We consider a longitudinal acoustic wave, which prop- 
agates along the SL i - axis (Fig. 1 inset), creat- 
ing a position and time (t) dependent potential energy 
field, V(x,t) = — U sin {k$x — cost), for each miniband 
electro n 32 ' 33 . The wave amplitude, U = eD, depends on 
the maximum strain, e < 0.5%, that the acoustic wave 
creates and on the deformation potential, D-. We con- 
sider acoustic waves whose wavenumber, ks, lies within 
the inner half of the minizone, so that there is linear fre- 



quency dispersion ws = v$ks, where v$ is the speed 
of sound. Since the sound wave exerts force along x 
only, the electron dynamics can be described by a one- 
dimensional model for motion in the lowest miniband. 
Within the tight-binding approximation, the kinetic en- 
ergy versus crystal momentum dispersion relation for this 
miniband is E(p. x ) = A[l — cos(p x d/h)]/2, where A is the 
miniband width, and d is the SL period^ 4 -. We take A = 7 
meV, d = 12.5 nm, D = 10 eV, and v s = 5000 m s _1 , 
corresponding to a GaAs/(AlGa)As SL used in recent 
experiment a 20 ' 34 , but obtain similar results for a wide 
range of SL parameters. The wave is sufficiently weak 
and spatially slowly varying to preserve the miniband 14 , 
thus ensuring the validity of a semiclassical model 3 ^. 
The semiclassical equations of electron motion are 



dx 



dt 



dH 

dp x 



Ad 



= = ^ = =sin ^ , (1) 
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dp x 
dt 



dH 

dx 



k s Ucos(k s {x + x ) - co s t), (2) 



where the Hamiltonian H(x,p x ) — E(p x ) + V(x,t). We 
solve Eqs. (1,2) numerically, taking v x = 0, and p x =0 
when t — 0, to determine the electron trajectories in the 
absence of scattering. 



A. Electron dynamics for initial position xq = 

In order to understand the general dynamics of a sin- 
gle electron, first we consider the simplest situation by 
setting xq — x(t = 0) = 0. 

We use the Esaki-Tsu modet^ 3 ^ to find the electron 
drift velocity 



Vd = (v x (t)exp(-t/T))/n 



(3) 



where (.) denotes integration over t > 0, taking, from 
experiment 2 ^, an electron scattering time r = 250 fs, 
which includes both elastic (interface roughness) and in- 
elastic (phonon) scattering. 

The solid curve in Fig. Q] shows Vd calculated as a 
function of U (lower scale), or, equivalently, e (upper 
scale) for an acoustic wave with u>s = 4 x 10 11 rad s _1 
and wavelength As = 2ir/ks ~ Qd less than the length 
of most SL a 17 ' 18 ' 40 ' 41 . For comparison, the dashed curve 
shows the usual Esaki-Tsu drift velocity^ calculated for 
an electron accelerated by a static electric field of mag- 
nitude kgU/e, where e is the magnitude of the electronic 
charge. As discussed extensively in the literatur e) 1 ' 2 ' 14 , 
the Esaki-Tsu Vd(U) curve is linear for small U, attains 
a maximum when U = Vet — h/r ~ 2.4 meV, and 
thereafter decreases with increasing U as more electrons 
complete Bloch oscillations before scattering. Figure [1] 
reveals that for both low and high U, the Vd(U) curves 
produced by the static force and acoustic wave converge. 
But for intermediate U there are major differences in the 
two curves. In particular, the acoustic wave generates a 
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FIG. 2: (a) Solid curve: electron trajectory, x(t), calculated 
for U = 4 meV. Within the white [gray] regions, V(x,t) is > 
[< 0]. Dashed line has gradient vs- (b) E'(p x ): dotted lines 
mark ±U when U = 4 meV. At arrowed peak, E (p x ) ~ Uc- 
Inset: V(x'), where dotted lines mark turning points of orbit 
in (a), (c) As (a) except U = 15 meV. Bloch oscillation bursts, 
within numbered brackets, are separated by sudden jumps 
(arrowed), (d) As (b), except U — 15 meV. Left- [right-] 
hand filled circles mark where E'(p x ) = U [-17]. Numbers 
label different minizones. Inset: adjacent wells (1 and 2) in 
V(x'), with arrows discussed in text. 

larger peak Vd value and a far steeper (factor « 13) NDV 
region. 

To explain these differences, we consider the electron 
dynamics in the absence of scattering. Figure Ufa) shows 
the x(t) trajectory obtained numerically from Eqs. (1,2), 
taking U = 4 meV, below the peak in the Vd(U) curve 
(solid curve in Fig. [I} generated by the acoustic wave. 
The trajectory consists of regular, almost sinusoidal, os- 
cillations superimposed on a linear background of gra- 
dient vs [dashed line in Fig. HJa)], suggesting that the 
acoustic wave drags the electron through the SL 42 i 43 . We 
confirm this picture by considering electron motion in 
the rest frame of the acoustic wave, in which the elec- 
tron's position, x'(t) — x(t) — vst, determines the static 
potential energy V(x') = — U sm(ksx'). In this frame, 
the Hamiltonian is H'(x',p x ) — E'(p x ) + V(x'), where 
E'{p x ) = E(px) — vsPx, and the equations of motion are 

, dx' 8H 1 Ad fp x d\ 

v * = Vx - Vs = lF = dpZ = 2H sm [-T)- vs > (4) 

dp x OH' 

~~dt = ^~dx F = cos ( k sx )• (5) 

Since H 1 is not an explicit function of t, it is a constant 
of the motion but does not equal the total energy, H. For 
the initial conditions considered here, H' = 0, meaning 
that E'(p x ) = —V(x') can only take values between ±Z7, 
marked by the horizontal dotted lines in Fig. 2(b) for the 



trajectory in Fig. 2(a). The lines reveal that the electron 
can only access the almost parabolic region of the E'(p x ) 
curve [thick in Fig. 2(b)] around p x = 0. Since for the 
given parameters the minimum value of E'(p x ) that the 
electron can attain is close to zero, its maximum potential 
energy is also close to zero. The electron is therefore con- 
fined within a single potential well in the acoustic wave 
and oscillates back and forth across this well between 
turning points at x' = and As/2 [vertical dashed lines 
in Fig. 2(b) inset]. Since the electron remains within 
the almost parabolic region of E'(p x ), where its effective 
mass is constant, x'(t) is an almost harmonic function of 
t. Therefore we can approximate x(t) as 

x(t) w«5t + As[l-cos(wflt)]/4, (6) 

where ljr is the frequency for motion to and fro across the 
potential well. This approximation accurately describes 
electron trajectories for small U , for example that shown 
in Fig. 2(a). The electron is trapped in the well, where 
V(x,t) < [gray bands in Fig. 2(a)]. But as the well 
moves, it drags the electron through the SL with a mean 
speed equal to v$ in the absence of scattering. 

Increasing U above 4 meV initially has no qualitative 
effect on the electron orbits. They continue to be dragged 
through the SL and are of the form x(t) = v$t + f(t), 
where the periodic function, f(t), becomes less harmonic 
as U [upper dotted line in Fig. 2(b)] increases, thus mak- 
ing the electron access nonparabolic regions of E'(p x ). 

When U reaches a critical value, U c , equal to the local 
maximum of E'(p x ) marked by the arrow in Fig. 2(b), 
the electron trajectories can reach the edge of the first 
minizone, and therefore change abruptly from closed to 
open orbits that can traverse several minizones. The lo- 
cal maximum of E'(p x ) occurs when dE(p x )/dp x = vs, 
i.e. when sm(p x d/h) = 2hvs/Ad (w 0.08 for the given 
parameters). Using small-angle approximations, it fol- 
lows that the local maximum occurs when p x w (far/d) — 
(2h 2 v s /A d 2 ) = p™. Therefore, from E'(p™) = U c , we 
estimate 

U c « A - v s hn/d. (7) 

Figure[5{c) shows x(t) calculated for U = 15 meV > U c 
= 6.2 meV. The bursts of high-frequency fluctuations in 
x(t) (within brackets) are Bloch oscillations driven by the 
acoustic wave. The jumps in x(t) (arrowed) occur at the 
centers of the white and gray stripes in Fig. 2(c), when V 
is extremal and, consequently, the acoustic force is zero, 
and therefore unable to induce Bloch oscillations. 

To explain fully the form of the trajectory in Fig. 2(c), 
we consider the electron motion in the rest frame of the 
acoustic wave. Initially, the electron is at x' = where 
the high gradient of V(x') [Fig. [2^d) inset] causes p x 
rapidly to increase to the edge of the first minizone [la- 
beled in Fig. 2(d)], thus reversing v x and v' x . After 
crossing the minizone boundary, the electron continues 
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to experience a large positive force, which increases p x 
through minizones 1-9 in Fig. 2(d), thus generating the 
Bloch oscillations within Bracket 1 in Fig. 2(c). As p x in- 
creases, the average value of E'{p x ) decreases [Fig. 2(d)] 
and V(x') increases (to keep H' — 0) as the electron 
moves up the left-hand side of Well 1 in Fig. EJd) in- 
set. As the electron climbs the well wall, \dV(x')/dx'\ 
decreases, thus reducing the frequency of the Bloch os- 
cillations and increasing their amplitudei^^, as shown 
by the x(t) curve within Bracket 1 in Fig. 2(c). 

When the electron reaches the top of Well 1, so that 
V(x') — U, E'(p x ) attains its lowest possible value of — U 
[lower dotted curve in Fig. 2(d)] and so p x can no longer 
increase. Instead, since the acoustic force is instanta- 
neously zero, p x is temporarily pinned at the intersection 
[right-hand filled circle in Fig. 2(d)] between E'(p x ) and 
the lower dotted line. The large negative velocity at this 
intersection, dE' /dp x f=s —5.6 x 10 4 m s _1 , makes the 
electron jump backwards along the section of the x(t) 
curve marked by the left-hand arrow in Fig. 2(c). This 
jump transfers the electron to the position marked by 
the right-hand arrow in Well 2 [Fig. [3](d) inset]. At this 
position, the acoustic wave exerts a large negative force 
on the electron, which causes p x to decrease, so inducing 
another burst of Bloch oscillations [within Bracket 2 in 
Fig. 2(c)], until E'(p x ) reaches its maximum value [upper 
dotted line in Fig. 2(d)] and V(x') attains its minimum 
value of — U in Well 2. Then, the electron again jumps 
backwards, along the x(t) trajectory marked by the cen- 
tral arrow in Fig. 2(c), with velocity sw —6.8 x 10 4 m 
s^ 1 , approximately equal to dE' jdp x at the intersection 
[left-hand filled circle in Fig. 2(d)] between E'{p x ) and 
the upper dotted line. This jump transfers the electron 
to the position marked by the left-hand arrow in Well 2 
[Fig. [2fd) inset], where a large positive force causes p x 
rapidly to increase, triggering the Bloch oscillation burst 
within Bracket 3 in Fig. 2(c). Thereafter, the cycle re- 
peats, with the electron jumping backwards after each 
Bloch oscillation burst. 

The number of Bloch oscillations within each burst 
equals the number of distinct minizones, N « 2Ud/hvs, 
that the electron traverses. When U = 15 meV, N = 14, 
corresponding to crossing the minizones labeled -4 to 9 
in Fig. 2(d). The abrupt onset of the acoustically-driven 
Bloch oscillations contrasts with the gradual switch on 
produced by increasing a static force 1 ^* 1 ^. 

Figure [3] shows a color map of the Fourier power, S(u), 
of v x (t) trajectories calculated for a range of U at fixed 
W5=4x 10 11 rad s _1 . The spectrum changes abruptly 
at U = U c (arrowed), due to the transition from the 
wave-dragging to Bloch oscillation regimes. 

For U < U c , S(u>) has a sharp peak (lower left light 
curve in Fig. [3j at the frequency, wr, for motion across 
the potential well [Fig. 2(b) inset] that traps the elec- 
tron and drags it through the SL. Three higher harmon- 
ics are also visible in the color map, but their power is 
orders of magnitude lower than the fundamental. When 
U ~ 4 meV, ujr » 17ws, indicating that the dragged 
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FIG. 3: (Color online) Fourier power spectrum, S(uj), (scale 
right) of v x {t) trajectories calculated for a range of U (lower 
scale) or s (upper scale) at fixed ljs — 4 x 10 11 rad s _1 . Ar- 
row is at U — Uc- Open circles [filled squares]: analytical 
estimates of ujr [ui^ ax ] in the wave-dragging [Bloch oscilla- 
tion] regimes. Solid and dashed white lines are discussed in 
text. 



electron paths cause significant frequency up-conversion 
of the acoustic wave. In the regime U < U c , correspond- 
ing to periodic x'(t) trajectories, the equations of motion 
yield a simple equation for ujr ps quj, where the factor 
a = (UA/n) 1 / 2 (d/hvs) is estimated by substitution of 
the approximation ([6]) into the set of equations (|4]) and 
([5]). This factor can be used to predict the frequency 
up-conversion attainable from a given SL. For the SL 
considered here, lor values obtained from the equation 
[open circles in Fig. [3J agree well with the numerically- 
calculated spectrum. 

When U exceeds U c , the bandwidth of S(oj) increases 
and the peaks become denser. The spectrum is broad be- 
cause the Bloch frequency changes continuously through- 
out each burst. The strongest peaks occur near the 
maximum frequency of the Bloch oscillations, u>g ax = 
ksUd/H^^, whose values are marked by the squares 
in Fig. [3J The series of abrupt jumps [arrowed in Fig. 
2(c)] between Bloch oscillation bursts generates the low- 
frequency (uj < 0.3 x 10 13 rad s _1 ) peaks in S(ui) and 
their harmonics. These peaks shift abruptly to lower ui 
as U increases (most easily seen for ui < 10 13 rad s~ 4 ). 
Two distinct series of jumps, each with a period of » 1.7 
meV, occur at U values marked by the solid and dashed 
white lines in Fig. [3] Their origin can be understood by 
considering Fig. 2(d). As U increases, the upper dotted 
line moves upwards through the E'(j> x ) curve. At the 
U values marked by the solid white lines in Fig. [3J the 
upper dotted line in Fig. 2(d) passes above a local maxi- 
mum in E'(p x ). This enables the electron to enter a new 
minizone, so adding an additional Bloch oscillation to 
each burst [within brackets in Fig. 2(c)]. As a result, the 
repeat frequency of the bursts decreases abruptly, thus 
red-shifting the corresponding spectral peaks in S(oj). 
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FIG. 4: (Color online) Color map of Vd versus U (or e: top 
scale) and A. Dashed line: U c versus A. us — 4 x 10 11 rad 




FIG. 5: (Color online) Color map showing the electron drift 
velocity, Vd, (scale right) calculated as a function of xo and 
U. Dashed curve shows U c (xo) obtained from Eq. Q and 
the dot-dashed curve is U' c {xo) calculated using {TUJ. 



Similar shifts occur at U values marked by the dashed 
white lines in Fig. [31 when the lower dotted line in Fig. 
2(d) passes below a local minimum in E'(p x ). 

In the limit U — > 0, u/r oc U 1 ^ 2 —J- 0, which means 
that the electron scatters when x' ~ [Fig. 2(b) inset] 
and so experiences an almost constant force, ksU. Con- 
sequently, in Fig. 1, the va curve for the acoustic wave 
(solid) converges to the Esaki-Tsu curve (dashed) for an 
electron accelerated by a constant electric field, ksU/e. 
The two curves also converge when U 3> U c because the 
electron immediately experiences a large positive force, 
equal to the maximum gradient of V(x'), which creates 
Bloch oscillations localised within a distance A /F (<C d) 
of x' = 0. Since U c > Uet (Fig. 1), the Vd curve pro- 
duced by the acoustic wave overshoots that generated 
by a static force and so causes a far higher maximal 
NDV value, Dy. Our analysis predicts strong acoustic 
enhancement of the peak Vd value, v T d nax , for all SLs with 
lj r t r; 1 when U ~ U c . This ensures that v^ ax is close 
to the mean speed (~ avg) of an electron traversing one 
well in the acoustic wave [Fig. 2(b) inset], rather than 
the lower speed, v$, of the well itself. 

Figure [4] shows a color map of Vd calculated versus U 
(or e) and A. For U > U c (dashed line), Vd decreases 
abruptly due to the sudden onset of Bloch oscillations. 
Figure 0] reveals that the velocity overshoot and, hence, 
Dy, both increase with increasing A. When A = 20 
meV, Dy ~ 60 times higher than for a static force, sug- 
gesting that wide miniband, acoustically-driven, SLs will 
exhibit very high frequency electron dynamics. 



B. Electron dynamics for initial position xo 7^ 

Since the propagating acoustic wave produces a 
spatially-varying potential it is important to consider 



how the electron's initial position, xq, affects its subse- 
quent motion. 

Figure^ shows Vd calculated as a function of U and xq. 
Due to the spatial periodicity of the propagating wave, 
Vd(xo,U) is also a periodic function of xq with a period 
equal to Ag. 

When xq = 0, corresponding to the Vd(U) curve shown 
solid in Fig. [TJ the color map in Fig. [S] confirms that Vd 
initially increases with increasing U > 0. However, when 
U = U c (x = 0) (vertical dashed line in Fig. [T]), the 
transition from the wave dragging to the Bloch regime 
produces a sharp suppression of Vd (red to yellow regions 
in Fig. [5]). Figure [5] reveals similar behavior for — As/4 < 
xq < As/8. By including xq explicitly in the Hamiltonian 
for the system, we find that H'(xq) = — U sin (ksxty), 
meaning that E'(p x ) = — U sin(fcsxo) — V(x',Xo) can 
only take values between ±U — U sin (ksXo)- Previously, 
we showed that for the electron to Bloch oscillate it must 
attain enough kinetic energy to traverse the first local 
maximum in the E'(p x ) curve [marked by the arrow in 
Fig. E[b)], which occurs when E'(p x ) rs A — vshn/d [see 
Eq. 0]. Therefore, 



U c (x ) 



1 - sin(fcsxo) 



VgtV7r\ 



(8) 



Hence, for U < U c (xq), the electron is dragged by the 
acoustic wave, whereas for U > U c (xq) it is allowed to 
perform Bloch-like oscillations. The values of U c (xq) ob- 
tained from Eq.® are shown by the dashed curve in Fig. 
[5] This curve is in a good agreement with the sudden sup- 
pression of Vd in the color map, which results from the 
onset of Bloch oscillations. 



For A* 



< 



■'Co 



< 3Ac 



the transition from the 



wave dragging to the Bloch oscillation regime occurs 
beyond the experimentally attainable range of U and 
there is no associated suppression of Vd- Instead, for 
given xq between As/8 and 3As/8, Vd is almost inde- 
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pendent of U . However, for fixed U, increasing xq pro- 
duces a gradual decrease in v^. This can be understood 
by considering the range of possible values of E'(p x ), 
given by ±U — U sin (ks%o), which implies that increas- 
ing Xo from to As/4 decreases the maximum attain- 
able value of E'(p x ) [upper horizontal dotted line in Fig. 
[DJb)]. As a result, the electron can no longer access 
the high gradient regions of E'{p x ) where the magni- 
tude of v' x = dE' (p x ) I dp x is high. Further increasing 
Xo causes v' x to decrease until Xo ps A5/4 at which point 
E'(Px) = 0. Consequently, the electron cannot oscillate 
within a potential well of the acoustic wave and is simply 
dragged through the lattice at a constant speed vs, so 
that Vd ~ Vs. 

For xo > As/ 4, the initial force on the electron 



dV(t = 0) 
dx 



= k s U cos(k s x ), 



(9) 



becomes negative. The electron therefore initially moves 
in the negative p x direction where the gradient of E'(p x ), 
and hence also v' x , is negative [see Fig. HJb)]. Therefore, 
Vd becomes increasingly negative as the electron starts to 
access the high (negative) gradient regions of E'{p x ). 

Figure [5] shows that when xo « As/2, increasing U 
from initially reduces ivz, increasing its magnitude, \vd\- 
However, at a critical value of U = U' c (xq), close to 
U c (xo), the magnitude of Vd < decreases dramatically 
(color map changes abruptly from dark to light blue in 
Fig. [SJ. One might expect that this suppression of \vd\ 
would occur exactly at the transition from wave drag- 
ging to Bloch trajectories, as seen when |x | < As/4 
where the initial force on the electron [Eq. ©] is pos- 
itive. However, for xq ~ As/2, the electron is initially 
forced in the negative p x direction [see Eq. fljjjl]. There- 
fore, for U w U c (xq), the electron scatters before it can 
traverse the local maximum in E'(p x ) [arrowed in Fig. 
[U^b)]. Consequently, when xq « As/2, the transition to 
the Bloch regime at U = U c (xo) has no effect on Vd- In- 
stead the suppression of \vd\ occurs when U is slightly 
larger than U c (xq). Specifically, the supression of \vd\ 
occurs when the electron can Bragg reflect by traversing 
the local maximum to the left of the origin at p x — 0, i.e. 
when E'(p x ) = A + vshw/d. Using this condition in H', 
we find that 



U' c (x ) 



1 - sin(fcsa;o) 



A- 



vshir\ 



(10) 



The values of U' c (xq) obtained from Eq. (TTU)) , shown 
by the dot-dashed curve in Fig. [SJ coincide almost ex- 
actly with the dramatic supression of \vd\ observed when 
x o = As/2. More generally, when — As/4 < xo < As/4, 
the initial force is positive [see Eq. ©], and U c (xo) accu- 
rately estimates the U value at which \i>d\ is suppressed. 
However, when As/4 < xq < As/2 and —As/2 < xo < 
— As/4, the initial force is negative [see Eq. ©], and 
U' c {xq) gives a better estimate of the position of \vj\ su- 
pression. Note that when Xq ss As/4, so that the initial 



force is and \vd\ is minimal, U c —U' c — > 00, meaning that 
the electrons never perform Bloch oscillations. 



III. COLLECTIVE ELECTRON DYNAMICS 
AND CHARGE DOMAINS 

To investigate how the acoustic wave affects 
the collective behavior of the electrons, we solved 
the current-continuity and Poisson equations self- 
consistently throughout the device. To do this, we 
adapted the widely-used drift diffusion model of mini- 
band transport in SLs— i 17 ' 39 , for the case of acoustic 
wave driving. In this model, we discretise the SL region 
into N = 480 layers, each of width Ax = L/N = 0.24 
nm, small enough to approximate a continuum. The vol- 
ume electron density in the m th layer (with right-hand 
edge at x = mAx) is n m and the electric field, F, values 
at the left- and right-hand edges of this layer are F m and 
F m+ \ respectively. In the emitter and collector ohmic 
contacts, F — Fq. The evolution of the charge density in 
each layer is given by the current continuity equation 



eAir- 



dn r , 



dt 



J m — 1 J n 



1 . 



.N, 



(11) 



where the areal current density from the m th to the m - 
1 th layer is 



en m u™ - eD 



dn m 
dx 



m 



1...N. 



(12) 



In Eq. (|12l) . the drift velocity in the m th layer, u™, is 
determined from Eq. ([3]), using the semiclassical equa- 
tions of motion (flj and ([2]), in which xq — mAx is 
the initial position of the electron, and the Hamilto- 
nian H{x,p x ) — E{p x ) + V{x,t) — eF m x where F m is 
the mean electric field in the m th layer. The diffusion 
coefficient, De, is calculated from the Einstein relation 
De — (fcflr/e)/i, where fi is the electron mobility in the 
linear part of the Vd{U) curve (see Fig. [TJi2, and the tem- 
perature T — 4.2 K. Since J m depends on the local drift 
velocity, v™(V, F m ), the collective electron dynamics de- 
pend directly on the single electron orbits. Conversely, 
the single-particle electron trajectories depend explicitly 
on the collective electron dynamics through their effect 
on the charge distribution and, hence, the local electric 
field F m . Our calculations include this mutual depen- 
dence because they are based on self-consistent calcula- 
tions of F m and of the single-particle electron trajecto- 
ries. 

In each layer, F m obeys the discretised Poisson equa- 
tion 



eAx 

F m +i = \n r 

€Q€ r 



n D )+F m , m = l...N, (13) 



where eo and e r = 12.5 are, respectively, the absolute 
and relative permittivities and tie — 3 x 10 22 m -3 is the 
n-type doping density in the SL layers^. 
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FIG. 6: (a) Upper (lower) solid curve: I ma x (/mm) calculated 
as a function of U. Dashed curve: time-averaged current, 

lave- (t*) la — Imax Imin VCXSUS U. 



We use ohmic boundary conditions, taken from a previ- 
ous experimen t ' 1 -^, to determine the current, Jo = <tFq, 
in the heavily-doped emitter of electrical conductivity 
a = 3788 S m _1 . The voltage (=0) across the device 
is a global constraint, which requires that = Vc + 
Em=i( f ™ + f m+i)j where the voltage, Vc, dropped 
across the contacts includes the effect of charge accumu- 
lation and depletion in the emitter and collector regions 
and a 17 f2 contact resistance^. Then, the total current 
flowing through the SL layers is I(t) = S m =o 
where A = 5 x 10 -10 m 2 is the cross-sectional area of the 

^14.17.39 

I(t) oscillates between minimum and maximum val- 
ues, / m i n and /max respectively, whose variation with U 
is shown by the upper and lower solid curves in Fig. [6^a) . 
The dashed curve in this figure shows the time-averaged 
current, J me . Remarkably, even at very low U values the 
current oscillates. This contrasts with the behavior of 
SLs with no acoustic wave driving, which exhibit current 
oscillations only when an applied bias voltage is large 
enough to induce NDV. Figure [U[a) reveals that as U 
increases from to 4 meV, the magnitudes of I max and 
Imin initially increase, with I max attaining a maximum 
value of pa 7.4 mA when U pa 4 meV, and I m i n hav- 
ing a minimum value of pa —5.8 mA when U pa 2 meV. 
The magnitudes of both I max and I m in decrease as U 
increases beyond ps 4 meV. 

In Fig. [BJb), we show how the amplitude of current 
oscillations, I a — I max — I m in, changes with U. Initially, 
I a increases with increasing U until it reaches a maximum 
of pa 11.6 mA when U ~ 3 meV. Thereafter, I a decreases 
with increasing U. The I(U) characteristics shown in Fig. 
IH1 can be understood within a single-electron picture. As 
discussed above, and shown in Fig. [5j all electrons follow 
dragged orbits when U < U c {xq — —As/4) pa 3.1 meV. 
Consequently, increasing U within this regime raises va 
(see Fig. [5]), thereby also increasing the magnitude of 
the current oscillations, since I oc Vd- As U increases 
beyond U c , Bloch oscillations gradually turn on, initially 
for electron trajectories starting at the maximum of the 
acoustic wave potential (with xq — —A/4) but eventually 




(ps) 



FIG. 7: I(t) curves calculated for U = (a) 1 meV, (b) 4 meV, 
(c) 15 meV. 



for the majority of electrons, thus suppressing the current 
oscillations. 

Figure [7] shows I[t) oscillations calculated for three 
acoustic wave amplitudes: (a) U = 1 meV, which corre- 
sponds to the ascending part of the I a (U) curve in Fig. 
|H{b); (b) U — 4 meV pa U c , which corresponds to maxi- 
mal I a in Fig. [Hb); (c) U = 15 meV, when I a is small. 
Remarkably, the frequency of the I(t) oscillations is in- 
dependent of U and equals the frequency of the acoustic 
wave uj s = 4 x 10 11 rad s^ps 63.7 GHz). 

Figure [7^a) shows that for small U, the I(t) oscilla- 
tions are almost sinusoidal indicating a single dominant 
harmonic component. However, as U increases the an- 
harmonicity of the I(t) oscillations also increases [see Fig. 
[7{b)]. For large U, above U c , the current oscillations are 
strongly anharmonic [Fig. [7{c)] due to the appearance of 
kinks in the I(t) profile. 

To understand the shape of the I(t) oscillations for dif- 
ferent U, we now examine the spatio-temporal evolution 
of the electron density, n(x,t), in the SL. Figure [8^ a) 
shows n(t,x) calculated for U = 1 meV within the wave 
dragging regime. The acoustic wave produces charge den- 
sity waves, which propagate through the SL with an ap- 
proximately constant speed (pa vs) and spatial profile. 
These charge density waves are a linear response of the 
SL to perturbation by the plane acoustic wave. In Fig- 
ure [5Jb), the solid curve (scale on left-hand axis) shows 
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f(ps) x/L 

FIG. 8: (Color online) Color maps of n(t, x) (left-hand col- 
umn) and instantaneous electron density profiles, n(x,t — 35 
ps), (solid curves in right-hand column) calculated for U = 
(a,b) 1 meV, (c,d) 4 meV, (e,f) 15 meV. Scale bars for color 
map are in units of 10 22 m~ 3 . Vertical dotted lines in (a),(c) 
and (e) indicate time t = 35 ps corresponding to plots in 
right-hand column. In (b),(d) and (f), dashed curves show 
the acoustic wave profile, V(x), at t = 35 ps (scale on right- 
hand axes) and solid (dashed) arrows mark electron density 
maxima (acoustic potential well minima). 



the spatial form of the charge density wave, n(t = 35 ps, 
x), i.e. along the vertical dotted line in Fig. IDJa). The 
dashed curve in Fig. |5fb) shows the potential energy 
profile, V(x,t — 35 ps), of the acoustic wave. Minima in 
the acoustic wave energy (dashed arrows in Fig. E^b)) 
lead to the local accumulation of electrons whose density 
is maximal [solid arrows in Fig. Efb)] near the acoustic 
wave minima and minimal near the acoustic wave max- 
ima. We find that the electron accumulation regions lag 
slightly behind the local minima in V: this is due to in- 
ertia as the electrons 'ride' up the left-hand sides of the 
potential wells as the acoustic wave propagates through 
the lattice from left to right. When each electron accu- 
mulation region reaches the collector contact [x = L) it 
produces a sharp increase in I(t). Another period of the 
charge density wave then forms near the emitter contact 
and the propagation process repeats, so producing I(t) 
oscillation s 14 ' 39 . For small U, the linear response of the 
electron gas to the acoustic wave means that the electron 
accumulation and depletion regions have similar spatial 
forms and magnitudes, making the I(t) oscillations al- 
most symmetrical around 7 = [Fig. [7[a)]. 



In the wave dragging regime, increasing U to 4 meV, 
increases both Vd and I a [see Fig. Efb)]. However, when 
the wave dragging force is combined with the electric 
field associated with the charge density modulation, elec- 
trons are occasionally driven into the Bloch oscillation 
regime, which localises them spatially. This localisation 
induces additional electron accumulation regions [light 
gray (cyan online) areas in Fig. [HJc)] , which are known 
as charge domains and are a nonlinear response of the 
electron gas to the driving forces. The charge domains 
appear as plateaux-like features in the (solid) n(t = 35 
ps, x) curve in Fig. HJd). This curve is less sinusoidal, 
and has higher peak values, than for [7 = 1 meV [Fig. 
IHKb)]. Consequently, the I(t) oscillations for U = 4 meV 
[Fig. E£b)] are both stronger and more anharmonic than 
when U = 1 meV [Fig. E£a)]. 

Further increasing U to 15 meV generally decreases Vd 
since the electrons are regularly driven into the Bloch os- 
cillation regime. The onset of Bloch oscillations causes 
electrons to localise and accumulate in certain regions 
of the SL, thereby creating the high electron density 
domains shown by the narrow dark gray (dark red on- 
line) stripes in Fig. [He). The associated electric field 
produces additional isolated islands of electron accumu- 
lation, shown by the light gray (cyan online) areas in 
Fig. HJe). In the (solid) n(t = 35 ps, x) curve shown 
in Fig. [8{f), these additional accumulation regions ap- 
pear as small broad peaks, which separate the sharper, 
dominant, maxima. 

Due to the co-existence of the distinct electron accu- 
mulation regions in different parts of the SL when [7=15 
meV [i.e. the small and large peaks in Fig. |5ff)], the I(t) 
oscillations [Fig. E£c)] have a complex anharmonic form 
with pronounced kinks just above I = 1.5 mA. Since 
the electrons are more strongly localised when U = 15 
meV, than when U = 4 meV, far fewer electrons per unit 
time arrive at the collector. Consequently, I max and I a 
are both smaller at the higher U value [see Figs. [HI and 
EP>,c)]. 

Figures |8ja),(c), and (e) reveal that a new charge ac- 
cumulation front forms shortly after each acoustic wave 
minimum arrives at the left-hand edge of the SL (x = 0), 
which occurs with a frequency us. Since this charge sub- 
sequently travels through the SL, 7(t) oscillates periodi- 
cally at a frequency ujs independent of the value of U. 

Note that the electric field resulting from electron re- 
distribution within the SL perturbs the single-electron 
trajectories, but has no qualitative effect on the crossover 
between the wave dragging [Fig. Ufa)] and Bloch oscil- 
lation [Fig. EJc)] regimes. However, the transition be- 
tween these regimes [Fig. [SJb)] is blurred slightly by the 
local electric field. This is because, near the crossover, in 
parts of the SL the local electric field opposes the electron 
acceleration due to the acoustic wave, thus maintaining 
wave dragging. But in other regions, the electric field 
supplements acceleration by the acoustic wave, so driving 
the electrons further into the Bloch oscillation regime. 
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IV. CONCLUSION 

In conclusion, acoustic waves can induce an abrupt 
transition between two distinct dynamical regimes of 
electron transport in SLs. When U < U c , the electrons 
oscillate within a single spatial period of the acoustic 
wave, at frequencies w# 3> wg, and are dragged through 
the SL with a drift velocity whose peak value can greatly 
overshoot that produced by a static field. For U > U c , 
the acoustic wave triggers bursts of Bloch oscillations, 
thus causing very high NDV. 

The acoustic wave causes the formation of charge den- 
sity waves and, when U > U c , nonlinear charge do- 
mains, which propagate through the SL and thereby cre- 
ate self-sustained current oscillations. The fundamental 
frequency of these oscillations equals that of the acous- 
tic wave for all U. However, both the shape of the I(t) 
curve, and the spatial profile of the charge domains, de- 
pend on U in a way that reflects the underlying single- 
particle dynamics. In the wave-dragging single-particle 
regime (U < U c ), increasing U strengthens the electron 
accumulation and depletion in the charge density waves, 
and hence increases the magnitudes of I max , Imin an d I a - 
Conversely, when U > U c , the onset of Bloch oscillations 
creates additional (nonlinear) charge domains, which co- 
exist with the (linear) charge density waves shaped by 
the acoustic wave period. In this regime, as U increases 



the electrons become increasingly localised, reducing the 
magnitudes of I max , Imin an d I a - 

The complex single and collective electron dynamics 
that we have identified create new perspectives for using 
acoustic waves to generate high-frequency electric cur- 
rent oscillations. They also highlight the potential of SLs 
to bridge the interface between conventional electronics 
and the rapidly-developing field of phononic a 28 ' 29 . For 
example, SLs could be used to transform acoustic waves 
into electro-magnetic ones, produced by the formation 
and propagation of charge domains, and, conversely, to 
detect propagating acoustic waves via the current oscil- 
lations that they induce. 

Finally, we note that cold atoms in optical lattices may 
exhibit similar dynamic a 30 ' 44 , with the abrupt transition 
to Bloch oscillations providing sensitive control of trans- 
port. 
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